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ABSTRACT 



The problem considered is the application of a one-dimensional 
maximum principle to second order, linear differential equations of 
the form 

u" + g(x)u* + h(x)u = f(x) for a ^ x ^ b 
with associated general boundary conditions to obtain functions z^(x) 
and z 2 (x) such that 

z 2 (x) ^ u (x) ^ z^x) 

on [a,b]. The functions f,g and h are assumed to be bounded. We wish 
to determine the behavior of the solution u(x) on [a,b] and also to 
obtain reliable numerical estimates of u. 

The basic concepts in the theoretical background are expanded 
versions of a presentation in Protter and Weinberger [Ref. 4]. 
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I. ONE -DIMENS TONAL MAXIMUM PRINCIPLES 



If u is a real-valued continuous mapping of a compact set X, then 
u(X) is compact and by the Heine-Borel Theorem, u(X) is closed and bounded. 
Thus, u is bounded. Then there exist points p and q in X such that u(q) ^ 
u(x) £ u(p) for all x in X; that is, u attains its maximum (at p) and 
its minimum (at q) . This is a specific example of a general maximum 
principle. We begin our discussion with a study of some specific maxi- 
mum principles; namely, where and under what conditions functions of a 
single real variable which satisfy differential inequalities can attain 
their maximum value. The approach parallels a similar investigation in 
Protter and Weinberger [Ref. 4]. 

2 

Definition : For u = (x) in C [a,b], define the linear operator L 

on (a,b) as follows: 

L[u] E U + 



where g(x) is any bounded function. 

2 

Lemma 1 : If u = u(x) in C [a,b] satisfies the differential in- 

equality L [u] ) 0 on (a,b), then u cannot attain a maximum at a point 
interior to (a,b). 

Proof: Assume, to the contrary, that u assumes a relative maximum 

at c, a( c ( b. Then u f (c) - 0 and u"(c) $ 0, hence 

LM = U?c)+ q(c)v!(c) = -u n (c)$0 



which is a contradiction of L [u ] y 0 on (a,b). Hence u can attain a 
maximum only at the end points a or b. We note that the boundedness of 
g was necessary in the proof of Lemma 1 since, if lim g(x) = ^ , then 
L [ u ] I would be undefined. The boundedness condition for g may be 
weakened. The fact that g(x) is bounded on every closed subinterval of 
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(a,b) will be sufficient for our purposes. It is important to note that 
g bounded on every closed subinterval of (a,b) does not imply that g is 
bounded on [a,b]. Consider g(x) = 1 / (1—1) on (-1,1), and g(x) unbounded 
as x — ► 1. 

2 

Lemma 2: If v * v(x) in C [a,b] satisfies the differential in- 

equality L [ v ] ( 0 on (a,b), then v cannot attain a minimum at a point 
interior to (a,b). 

Proof: Note that -v satisfies the hypotheses of Lemma 1 since 

differential operators are linear; that is L [-v] = -L [v] > 0. 

Hence -v cannot attain a maximum interior to (a,b) which implies v cannot 
have a minimum interior to (a,b). We now have as tools two simple one- 
dimensional minimum and maximum principles. For convenience we will 
generally discuss explicitly only the maximum principles; but with 
reasoning similar to the proof of Lemma 2, we will be implicitly dis- 
cussing the related minimum principle. 

To continue our development we relax the condition of strict in- 
equality on the differential inequality and now consider functions u 
such that L [u] 0. The solution u identically a constant does not 

satisfy the strict inequality in the hypotheses of Lemmas 1 and 2 since 
L [u]HQ, We now allow as a special case the constant function, which 
attains its maximum at every point in [a,b]. We will show in a Corol- 
lary to Theorem 1 that the functions u H constant are the only func- 
tions which satisfy L [u]^0 and attain their maximum interior to (a, 
b). 

In investigating the behavior of a function u which attains its 
maximum at x = a, we might imagine that u could have a graph of the form 
shown in Fig. la (and similarly if the maximum is attained at x = b). 
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That is, u attains its maximum at x = a with u'(a) = 0 (or at x = b 
with u fl (b) = 0) . 






(a) (b) (c) 



Fig. 1 

We show now that this behavior is impossible if u satisfies L [u] ^ 0, 

If M is the maximum of u on [a,b], then u(a) = M implies u f (a) <C 0 

(Fig. lb): and u(b) = M implies u'(b) 0 (Fig. lc) . 

2 

Theorem 1 : Let u - u(x) in C [a,b] be a nonconstant function 

which satisfies the differential inequality 

V'+ >, o 



on (a,b) and has one-sided derivatives at a and b, and suppose further 
that g is bounded on every closed subinterval of (a,b). Then: 

(i) the maximum of u occurring at x = a and g bounded below at 
x - a imply u 1 (a) ^ 0; 

(ii) the maximum of u occurring at x = b and g bounded above at x ~ 
imply u’(b) ^ 0. 

Proof: Suppose that u(a) = M and u(x) ^ M for a < x ^ b. Since 

u is nonconstant by hypothesis, there exists a point d in (a,b) such 
that u(d) ^ M, Define the function 



ztn) ~ 



A 0C-*') 

e 



i 



b 



with A a positive constant to be determined. We note that z(x) is 
positive on (a»b) and that 
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Lte)= - A 2 e 



+ goO/Ae 



= A [ A + 9^)J 



a.) 

e 



We choose A large enough so that L [z] ) 0 on (a,d); that is, we 
select A such that it satisfies the inequality 

A > - 9<« 

on [a,d] which is possible by the boundedness assumptions on g(x). 

Now form the function 

wa*) * + E Zft 1 ) 

where E is a positive constant chosen so that it satisfies 

I < A -> 14 > 

*(4> 

This is possible since the numerator is positive and since u(d) ^ M. 
Now 

LM = LM + EU*]>o 

on (a,d) and hence the maximum of w on [a,d] must occur at one of the 
end points by Lemma 1. We have 

V(&) - U(a) 4- E ?(A) - H 

and 

W(4) « *(4) + £ < V<A) +JH- - J*L 

upon substituting the inequality for E; hence w(a) ^ w(d) implies 
the maximum occurs at a. Obviously, the one-sided derivative of w at 
x - a cannot be positive for then w(a) would not be a maximum. Thus 

* Tl'(a) 4- E ifa) 4 O 
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and since z 1 (a) = A } 0, it follows that u'(a) ^ 0, the desired re- 
sult in (i). 

Suppose that u(b) = M and u(x) ^ M for a x ^ b. Since u is 
nonconstant by hypothesis there exists a point d in (a,b) such that 
u(d) ^ M. Define 

ZOL) » & - 1 

with A a positive constant to be determined. Then z(x) is positive on 
(a,b) and 



L0] * 2' V gflo 2' - 



> -A flL-b) 
A e 



A gMe 






= A [ A - 0W 



- Afa- fc) 

e 



We choose A large enough so that L [z] ) 0 on (d,b); that is, we 
select A such that it satisfies the inequality 

A > 9«.) 

on [d,b] which is possible by the boundedness assumptions on g(x). 

Now we proceed as before and form the function w(x) = u(x) + E z(x),and 
we reach the conclusion that w attains a maximum at x = b . The 
one-sided derivative of w at x = b cannot be negative for then w(b) 
would not be a maximum. Then 

w'(b) = H'(b) + E i '&) )0 

and since z'(b) - -A < 0, it follows that u ! (b) > 0, the desired re- 
sult in (ii). 

Corollary iA : (One -Dimensional Maximum Principle) 

2 

Suppose u - u(x) in C fa,b] satisfies the differential inequality 
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u" + g(x) u 1 } 0 on (a f b), with g(x) bounded on every closed subinterval 
of (a,b). If u attains the relative maximum M at an interior point c 
of (a,b) f then u(x) 2 M. 

Proof: Assume u has a maximum at c in (a,b) f then u'(c) = 0. Ap- 

plying Theorem 1 to the interval (a,c) and (c,b), we conclude that u is 
constant. This concludes the proof. We may prove Corollary 1A in- 
dependently of Theorem 1 as follows: 

Assume that u(c) = M and that there exists a point d in (a,b) such 
that u(d) < M. Take d > c for convenience. Then define the function 

A A 

zm = (p-a.) - (c-a) 

where A is a positive constant to be determined. We note that: (i) 
z(x) ^ 0 for a < x < c, (ii) z(c) = 0, and (iii) z(x) > 0 for c < x 

< b. We have 

LB = z" + 2 ' 

4-2 A- 1 

* A(A-0(%-*) + gooAO*-«0 

P A-Z 

* a[(a-<) -4- 9(A) a-*)] it-*) 

We choose A large enough so that L [z] ^ 0 for a^ x ^ d; i.e. f we 
take A to satisfy 

(A-i) + (X-- a)> o 

or equivalently 

A > - $(*■)(*■-*) + 1 

which is possible since g is bounded. Define the function 

= MW +■ E 

where E is a constant satisfying 
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o < E < 



JA_- K(d) 



The existence of E is guaranteed by the hypothesis on the point d and 
the function z. Hence on (a,c) 

W(t) - Idtt) + JE 2 ft) < A 

by (i). Also 

t/j/A) s? y(d) + B 2/d ) </ y/d) / At ~ u/d) =- JdL * 

At the point c 

'HJ(C) ~ V(C) + E & Ji . 

Hence w attains a maximum greater than or equal to M at an interior 
point of (a,d). But 

LM * + E Lffl > ° 

by construction of the function z, and we contradict Lemma 1. If 
d < c, then use the function 

A A 

2(H) - (c- a) - ft-*) 

which satisfies (i) z(x) ) 0 on (a,c) (ii) z(c) = 0, and (iii) z(x) 

^ 0 °n (c,b). Note that the inequalities (i) and (iii) are reversed 
for this case. This concludes the proof. It would also be possible to 
use the function .... ^ 

AilC-c) 

2ft) = e - 1 

as the auxiliary function in the proof of Theorem 1, if d ^ c, and 

„ . -AU-c) . 

&*■>= e -l 

if d ^ c. 
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Corollary IB : A nonconstant u satisfying the differential in- 

equality L [u] ^ 0 on (a,b) cannot have a relative maximum at an 
interior point. 

Proof: Assume, to the contrary, that u is nonconstant and has a 

relative maximum at an interior point c. Then apply Corollary 1A to a 
subinterval I containing c, where I is small enough that u(c) is an 
absolute maximum on I. This yields the immediate contradiction that 
u 5 M. 

Corollary 1C : A nonconstant function u satisfying the differential 

inequality L [u] ^ 0 on (a,b) can have at most one relative minimum in 
the open interval (a,b). 

Proof: We note that between any minima there must be a relative 

maximum which contradicts Corollary IB. 

The boundedness properties of g are required for Theorem 1 and 
Corollary 1A as shown in the following example. Consider the functions 
u(x) * cos x and g(x) = -cot x, which satisfy the differential equation 

M" + ffrjv* t£ - CoS i - cd t (- *-) 

- —Oast + cozt - O 

and here u M + g(x) u 1 ^ 0. Conclusion (i) of Theorem 1 does not hold on 
[0, / TK /2) since u f (0) = -sin(0) = 0 and u is a nonconstant function. 

On [- IK/2, IK/2) the conclusion of Corollary 1A fails to hold since 

u(0) = cos(0) = 1 is a maximum at an interior point and u is a noncon- 

stant function. We note that g(x) is unbounded in the neighborhood of 
the origin. 

Definition : A function u(x) has a horizontal point of inflection 

at x = c if u’(c) = 0 while u is strictly increasing or strictly de- 
creasing in some interval containing c. 
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Corollary ID : A function u satisfying the differential inequality 



L [u ] £ 0 on (a,b) can have no horizontal point on inflection in (a,b). 

Proof: Assume, to the contrary, that at the point c in (a,b), u f (c) * 

0. Then on [a,c] or [c,b], u must attain its maximum on the interval at 

x = c, a contradiction to Theorem 1. 

2 

Definition : For u = u(x) in C [ a ,b ] and h(x) and bounded function, 

define the linear operator (L+h) on (a,b) as follows 

CL + h)[u]E.'H* + h/jc)U 

where g and h are any bounded functions. 

In this more general setting, it is necessary to modify some of the 
preceeding discussion as shown by the following examples. Consider the 
differential equation 

u"+ ^ = 0 on 

which has the solution u(x) = sin x. Now u attains a maximum of 1 at 
X = 1 ^/ 2 , an interior point of [0,1l^]. 

We also see that the condition h(x) ^ 0 is not sufficient to yield 

r ' 

an unrestricted maximum principle. Consider 

Ti ~ O Cn [~/j /] 

X “X 

which has the solution u(x) = -e -e , and u attains a negative maximum 
of -2 at x = 0. We will prove that a nonconstant function u satisfying 
(L+h) [u] 0 with h(x) ^ 0 on (a,b) cannot attain a non-negative maxi- 

mum at an interior point. 

Lemma 3 : If u satisfies (L+h) [u] ^ 0, with h(x) ^ 0 on (a,b), 

then a nonconstant function u cannot attain a non-negative maximum at 
an interior point. 

Proof: Assume, to the contrary, that u(c) = M for a ^ c < b and 

u(x) ^ M on (a,b). Then h(c)u(c) ^*0, u'(c) = 0, u M (c) £ 0; hence 
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(L+ « + y(e) it'(c) + h(c)7t(c) 




- Tt'Yc) ■+ hfc) n(t) ~ ° 



a contradiction of (L+h) [u] ) 0 on (a,b). 

2 

Theorem 2 : If u = u(x) in C [a,b] satisfies the differential 



inequality 



(L+h)M ~ 7t"+ ${*■)*' v- ha)H ^ O 



on an interval (a,b) with h(x) ^ 0, if g and h are bounded on every 
closed subinterval, and if u attains a non-negative maximum value M at 
an interior point c, then u(x) 22 M. 

Proof: Suppose that u(c) = M and u(x) £ M for a < x < b and M £ 0. 

Then assume that there exists a point d in (a,b) such that u(d) M and 
for convenience take d ^ c. Define the function 



where A is a positive constant to be determined* We note that: (i) 
z(x) ^ 0 on (a,c): (ii) z(c) = 0; and (iii) z(x) ^ 0 on (c,b). Now 




(L+h) Li] * S>/*J ? 



* A e + AyMe + 
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We choose A large enough so that (L+h) [z] ^ 0 on (a,d), that is, 

so that A satisfies the inequality 



Since 



[a 1 * Agfa) 4 htti] - hit) 






>0 and h =: 0 on (a,b). 



-AU-c) 

e > o. 

we need only to choose A such that 



A$(t) + h(t) > o 

which is possible since g and h are bounded. Now form the function 

w; - im /- a i/f.) 



where E is a positive constant so that 

_ _ M- ftf'D 
e < *(<i) 

Hence on (a,c) 

TV/*) - *()tt + £ 2/X) < JL 



since z(x) ^ 0 on (a,c). And by the definition of E, 

TP/d) = u(d) -h M- 74/4) » AL 

Now 

Me) - t £ -2/c) - M 

since z(c) = 0; thus, we may conclude that w has a maximum greater than 
or equal to M which is attained interior to (a,d). But 

(Z*h)[w] = + Z (z+Wil > 0 

on (a,d) and we have a contradiction to Lemma 3 since w(x) is noncon- 
stant. 



If d ^ c, construct the function 

a/K) j 

*r£j ~ e. - 1 

which has the properties: (i) z(x) ^ 0 on (a,c), (ii) z(c) = 0, (iii) 

z(x) ^ 0 on (c,b). Now 
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(L +h)te 1 - £ " + h/xj 2 



-Ate-*) r-Ate-c) -> 

- Ay*)e +hte)[e - 1J 
- e Aj(t) + hh - hte) 








We choose A large enough so that (L+h) [z] ^ 0, that is so that the 
followine ineaualitv is satisfied: 



which is possible since g and h are bounded. We note that the condition 



will suffice for both d c and d ^ c. The remainder of the proof is 
exactly the same as before with c ^ d. 

Corollary 2A : If u satisfies (L+h) [u] ^ 0 on (a,b) with 

h ^ 0 and h(x) 0, and u assumes a non-negative maximum value at an 
interior point, then u(x) 0. 

Proof: Let H ^ 0 be the maximum attained by u on |a,b], By 

Theorem 2, u(x) £ M and hence 



since £ 




, choose A such that 



A - >o 



that A satisfies 



/j L - ft /$i*) / -I- k/x) > C 



(J>*h)L#Q - My 0 
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But h(x) £ 0 and h(x) 0 implies M = 0. 

2 

Theorem 3 : Suppose that u = u(x) in C [a,b] satisfies the dif- 

ferential inequality 

(L+h)[\ ]* K"+ + hWH ± O 

on (a,b) and has one-sided derivatives at a and b, that h ^ 0 and that 
g and h are bounded on every closed subinterval of (a,b). Then: 

(i) the non-negative maximum of u occurring at x ‘ a and the func- 
tion g(x) 4 - (x-a)h(x) bounded from below at x s a imply that u*(a) 0; 

(ii) the non-negative maximum of u occurring at x = b and the func- 
tion g(x) - (b-x)h(x) bounded from above at x 3 b imply that u' (b) ^ 0. 

Proof: We may extend the proof of Theorem 1 by noting that if the 

maximum of uc occurs at x = a, then 

(L+h)[z -j Q*e Ja + #$(*) + l- c )J 



± e 



am-*) r 



' 4 z y- A gw + A Oc-a) h/*)J 



i 

since Q £ TT" for t ^ 0 and h(x) Z 0 by assumption. If g(x) + 

•£+/ 

(x-a)h(x) is bounded from below at x = a, then we may choose A such 



that 



A > - [s (*■) + (*■-*) httf) 



If the maximum of u occurs at x = b, then as above 



(Mlz*-*- 1) 






- C 



’ [ft - A 3 (x) 1 iJ )J 






_ A/X-tJ 



-t A fb-x) 
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Hence if g(x) - (b-x)h(x) is bounded above at x = b, then we may choose 
A such that 

A > - x) hte) 

The remainder of the proof is analogous to Theorem 1. 

Corollary 3A ; If u ^ 0 satisfies (L4h) [u] ^ 0 on (a,b) with 
2 

h 0, if u is in C [a,b], and if u(a) & 0, u(b) ^ 0, then u(x) ^ 0 
on (a,b). 

Proof: Assume that there exists a point c on (a,b) such that u(c) 

^ 0. Then u(a) ^ 0 and u(b) ^ 0 imply that u attains a non-nega- 

tive maximum at some point d interior to (a,b),say u(d) = M. Then, by 
Theorem 2, u(x) S M ^ 0, a contradiction. 

A useful technique in determining the existence of a maximum (or 
minimum) of a solution to a differential equation at an interior point, 
independent of the preceeding theory, is a simple direct argument. 
Consider the problem 

a"* = — x on foiOj 

which does not have the format of the problems we have so far considered 

K 

due to the exponential term £ . 

We note that 

'f(y) s k y V / < o 

on (0,1). Hence, u cannot attain a minimum in (0,1); otherwise, if 
there exists a point c in (0,1) such that u(c) = M and u(x) ^ M for 
x ^ c; then u"(c) ^ 0 and we have 

on (0,1), a contradiction. 

As another example of the preceeding, consider 

( fy') =■ £ 
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whose solution u cannot attain a local maximum in the interior of any 
interval. Let I be an interval with subinterval [a,b]. If there exists 
c in (a,b) such that u attains a maximum in [a,b] at c, then 

U"(c) ±0 and K'/O-b-O 

Hence 

74 n fc) = i+ 2 cosfute)] - 2 Cos /t>) = 2 } 

an impossibility. 

We may even apply the previous results to boundary value problems. 
Consider 

n + c t< « -± 

for Oxl, with u(0) = u(l) = 0. Now 

e\' < o 

on (0,1); hence u cannot attain a minimum interior to (0,1) as in the 
first example. The boundary conditions u(0) = u(l) = 0 imply that u f (0) 
^ 0 and u f (l) 0; otherwise u would attain a minimum interior to 
(0,1). If u f (0) = 0 and u(0) = 0, then from the differential equation, 
u"(0) = -1. The boundary condition u(l) = 0 forces u to attain a mini- 
mum interior to (0,1), a contradiction. We obtain a similar result if 
we assume u(l) = u f (l) = 0, since u(0) = 0. Thus u* (0) 0 and u'(l) 

< o. 
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II. A GENERALIZED ONE -DIMENS TONAL MAXIMUM PRINCIPLE 



If we drop the restriction that the function h must be non-positive 
we may then obtain a maximum principle for (L+h) [u] ^ 0 by introduc- 
ing a function w, with certain desirable properties, such that the 

quotient u/w satisfies the previous maximum principles* 

2 

Theorem 4 : Suppose u = u(x) in C [a,b] satisfies the differential 

inequality (L+h) [u] ^ 0 on (a,b) with h(x) bounded and g(x) bounded 

below; then for any sufficiently small subinterval [a^b*], there exists 
2 

a function w in C [a,b] such that: 

(i) w > 0 on [a ' ,b 1 ] 

(ii) (L-fh) [w] ^ 0 on (a^b 1 )* 

Then the function u/w satisfies the maximum principles of Theorems 2 and 
3 on (a* ,b ') . 

Proof: Consider the function 

(M/t) * i - «') 

for a ^ a 1 ^ x ^ b* ^ b, where B is a positive constant to be 
determined. When B has been determined, we will suppose that (b^a 1 ) 
is such that 

(i) B(b’-a , ) 3 '< t 

so that w ) 0 on [a'jb 1 ]. Now 



(L 5 w"+ ymiv* + bit) uJ 

= (-2 it) + 

» -2r$j-i-h + { (*-*')* h/t)J + bit) 
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Since g and h are bounded below, choose constants G and H such that 
g ^ G and h ^ H. We now make a restriction on (b^a 1 ) so that 

(2) i * fo-ft') + { (*’'*') H > o 

for a' 4 x £ b' and choose B so that (L4h) [w] £ 0 on (a',b'), 

i,e * » ( hit) 

S " si u 

a permissible choice since h is bounded above ♦ Hence if the subinter- 
val [a^b 1 ] is sufficiently small, the function w with properties (i) 
and (ii) may always be constructed. 

Define the function v “ u/w on [ a f ,b * ] ; then 

(U+h)[Tij* f/Lr+h)[™] - + hlfunt 

= /yui"+ 2v'vs' i-u/Y ") VU /Vm w'it) 

Now dividing by the positive function w, we obtain 

r" + { 2 M + + r - °- 

Thus, v satisfies the differential inequality 

r"+j,r»v' a h,ft)v i o 

on (a'jb 1 ) where 

g, ft) = 2 £ + yr) ; ^ ft ) » . 

2 

By the construction of w, w in C [a'jb 1 ]; thence g and h are bounded 
on [a f ,b ? ]. Also h is non-positive by conditions (i) and (ii). Thus 
v satisfies the hypotheses for Theorems 2 and 3 on subinterval [a'jb 1 ], 
which concludes the proof. 
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We now see that a function u which satisfies (L+h) [u] ^ 0 has 



certain restrictions on its zeros on the interval under consideration. 
Corollary 4A : If the function u satisfies Theorem 4 on (a 1 ,!* 1 ) 



negative) on (a 1 ,!* 1 ). 

Proof: We may take the interval [a 1 ,!* 1 ] such that u(a*) = u(b*) - 



Then v(x) 55 u(x)/w(x) > 0 for a 1 x ^ b 1 and v(a ! ) = v(b') = 0. 

Hence v has a non-negative maximum in (a^b 1 ) which is a contradiction 
to Theorem 4. Thus, u(x) ^ 0 on (a',b'). Also, there could not be more 
than two zeros; otherwise, if there exists a point c in (a^b 1 ) such that 
u(c) * 0, then u(x) 0 on (a f ,b') forces u f (c) = 0. But then u must 
be the trivial solution, a contradiction. 

Note that the zeros of u must be isolated; i.e., if u(a) = 0, then 
u cannot vanish on some finite interval to the right of a. 

Corollary 4B : If u is a solution to the differential equation (L+h) 

[u] * 0, then u can have at most one zero in any interval (a'jb 1 ) where 
Theorem 4 holds. 

Proof: Note that -u also satisfies (L+h) [-u] =* 0. By Corollary 

4A, u and -u can each have at most two zeros on (a^b 1 ), between which 
both u and -u must both be negative, an impossibility. Thus, u has at 
most one zero on (a^b 1 ). 

We may conclude that the boundedness hypothesis in Theorem 4 on the 
function g is essential by noting the following example. The function 
u(x) “ x satisfies u M + g(x)u* = 0 with g(x) = -2/x, and g is unbounded 
in the neighborhood of zero. Consider the interval of the form (a,b) 
with b * 0. Then u attains a non-negative maximum of zero at x = b and 




then u can have at most two zeros (between which u is 




0 for a 1 ^ x ^ b. Assume that u(x) 0 for a* ^ x^ b. 
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u*(b) = 0, contradicting conclusion (ii) of Theorem 3, hence contra- 
dicting Theorem 4. 

Definition : Let r(x) be a solution of the differential equation 

r " gfc) r' v- h/u f - o 

with g and h bounded functions. If r 5 / 0 and r(a) = 0, the r must be 
non-zero on some finite interval to the right of a by Corollary 4B. If 
r has another zero on (a, ), denote the first such zero by a* and call 

it the conjugate point of a. (If r / 0 on (a, <=*>), then set a* = ) . 

Thus r has one sign on (a, a*) and for convenience we may take r(x) 

0 for a x^a*, since -r satisfies the same differential equation. 
We now establish a necessary and sufficient condition for the existence 
of the function w of Theorem 4 in terms of a*. 

Corolllary 4C : If a* is the conjugate point of a, there exists a 

function w 0 such that Theorem 4 holds on the interval [a,b] if and 
only in b a*. 

Proof: Assume that b ^ a*. If w ^ 0 on [a, a*], then the func- 
tion u/w vanishes at a and a* since r(a) * r(a*) = 0 , and u/w is posi- 
tive on (a, a*). Hence u/w attains a non-negative maximum in (a, a*) and 
u/w is non-constant, then w cannot satisfy (L4h) [w] ^ 0 on (a, a*) with- 
out contradicting Theorem 4. Hence by a contrapositive argument, we 
have the forward implication that the existence of an acceptable w(x) 
implies b ^ a*. 

To get the reverse implication, let b in (a, a*). Note that r(x) is 
bounded from below by a positive number m 3 m(c) on any subinterval [c,b] 
with a < c. Let 

s-m - 2 - e 

where A is a positive constant to be determined. Consider the function 
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m* j c wj + & zsjl) 

Now choose E small enough so that w ^ 0 on [a,b]. Now 

(L^}+ E(l 1 ih)[z-c A, *'* ) ] 

• e {(-&***> * yv(-*ryhM( 2 -c*n} 
■ - e. Mn, fs*» ,» > t,»[l - ti * i 



We want to choose A such that (L+h) [w] 6 0 on [a,b] which can be done 
as follows: 

AO- -a) A (Ad) 

Now -Ee ^ 0 since E and e are positive. We make 

two restrictions on A: 

-am-*) , /nZ 

(i) 1 - 2e >0, i.e., ft * ' 

Since h is bounded, we may also choose A large enough so that 

r _ f#~A) ~] 

ai) /j z -h Aqcn) ■+ hfx)ji-2£ j > #' 



Then we have (L+h) [w] 4s* 0, and w is a function for which Theorem 4 holds. 
If h(x) is unbounded or g(x) is not bounded below there may be no 
interval [a^b 1 ] for which Theorem 4 holds. Consider the function 

7i(t) * K 

which satisfies 

W xT%c ~ o 

Now as x — > 0+, h(x) = x — *■ + •**. There can be no function w ^ 0 
with the property that u/w satisfies the maximum principles in any of 
the intervals [0, l/n'iK ] (n = 1,2,...) since u(l/n^T') = 0 and l/nlr^O 
as n — ► 00 . Thus, for all £ > 0, the interval (0, £ ) contains a 

zero of u and Corollary 4A denies the existence of an appropriate w(x). 
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As in the examples in the previous section, we use a direct argument 



to conclude that no solution of 

// ^ 

# t e s 

can attain a minimum value greater than 1 or a maximum value less than 1. 
Assume there exists a point c such that: 

(i) u(c) = maximum; hence u n (c) 0 and we have 

rtte) 

e - e = - taHv) ^ 6 

KtCi 

which implies e > e, and we have the maximum value of u is u(c) ^ ln(e) 

- i; 

(ii) u(c) = minimum; hence u"(c) 0 and we have 

U (&) 

e 

UCC) 

which implies that e 
u(c) £ ln(e) = 1. 



e « u'Yc; ^ o 

e and we have the minimum value of u is 
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Ill . INITIAL VALUE PROBLEM 



Consider the solution to the differential equation 

(1) (t) ti'+ hum - -f 

which satisfies the initial conditions 

( 2 ) MV)*#, ) K'fa) = % 

where the functions f, g and h are defined on the interval (a,b) with 
g and h bounded and ^ and ^ given constants. 

The existence of a solution follows from the classical theory of 
ordinary differential equations [Ref. 3]. We may prove the uniqueness 
of such a solution independent of the classical theory using Theorem 4. 

Theorem 5 : Suppose the u^(x) and u^(x) are solutions of (1) which 

satisfy the same initial conditions (2) in an interval (a,b). Then 
n l = u 2 in (a,b). 

Proof: Define u(x) * u^(x) -u 2 (x). Then u satisfies 

7t‘ J + htfi) it *■ 0 

with initial conditions 

?t (a) * vW ~o 

Now assume that u ^ 0 in (a,b). Then by Theorem 4 there exists a 
positive function w on an interval (a, a + £ ) such that the maximum of 
u/w in (a,a+£ ) must occur at one of the endpoints. We note that 
is a positive constant depending only on the bounds of g and h. I.e., 

(a+ £ ) - £ * £ must satisfy conditions (1) and (2) in the proof of 
Theorem 4* But -u also satisfies the same equation with the same initial 
conditions, hence, also by Theorem 4, -u/w must also attain its maximum 
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at one of the endpoints, a or a+. Hence u/w must attain either a maxi- 



mum or a minimum at a. But at x = a 



(- 

V \A) 






U'W- MW 1 

u> z 



*=*0 



K 



since u(a) = u 1 (a) = 0. Now the function u/w satisfies Theorem 3; hence 
u/w must be identically a constant, and from the initial condition, u=0 
on [a,a+<2 ]. In particular 



Hence, by the analogous argument, u sEr 0 on [a+£ , a+2£. ] with £ un- 
changed since the choice of £ is independent of u. Now (b-a) < ; 

hence, after a finite number of steps, we may conclude that u = 0 on 

[a,b]. Thus u^S. u^* 

Consider the problem 

M" — U* = o 

with the initial conditions 

?/(o) - -H Vo) = O 

_ 2 
We note that u^ ~ 0 satisfies the equations trivially and for u^ =* x , 

i ' = 2 - f ~ 0 ■ 

But u 1 ^ on (0,1). Since g(x ) — >- as x — >~ 0, we may conclude 

that the boundedness of g is necessary in Theorem 5# 

We can easily find a function h(x) such that the equation 

% /; •+ H tf) & = o 



has two solutions satisfying 

7/ So) - # Ud) ~ O 



Obviously u ; S. 0 is a solution. We must then pick h(x) unbounded in 

the neighborhood of the origin; otherwise Theorem 5 would imply that the 

2 

trivial solution was the only solution. Choose h(x) = -6/x , then u(x) 

3 

= x is an appropriate solution satisfying the initial conditions. 



Hence we see that the boundedness hypothesis on h is also necessary in 
Theorem 5. 

An interesting application of the Mean Value Theorem arises in 
the initial value problem 

^ •=■ 4 Y X > O 



(3) ^ •+ e 



satisfying 



(4) 



7//6J = / and t dto) = 0 . 



Assume that and are solutions of the differential equation (3) 
satisfying the same initial conditions (4). Then the function 

7//XJ = % (X) - 74 /X) 

satisfies the homogeneous differential equation 
By use of the Mean Value Theorem 

tt, u? *% + eWi-Ui) 

e - e = f d,-u 2 ) e 

and the differential equation becomes 

7h+ OfUt-tii) 



; 0 <o^l 



with 



d(0) ~ 74' to) -o 



Hi * 6 (H, - Hi) 

Since h(x) = -e ~ 0 and is bounded on any finite sub- 

interval of (0, ), we may apply Theorem 5. The trivial solution 

u S' 0 satisfies the above initial value problem, thus u is unique 
and u^ sZ u^. Hence we have the solution of (3) satisfying (4) • 
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IV. BOUNDARY VALUE PROBLEM 



We will first be concerned with boundary value problems of the 

type 

(i) j (*■)%' +htt) w - £{%) 

on the interval (a,b) subject to the boundary conditions 

( 2 > iff#) ~ If, ') v.(io) ~ ' . 

We may derive a uniqueness theorem by the use of Theorem 4 for solutions 
of (1) which satisfy (2). The sign of the function h(x) has great signi- 
ficance as seen in the simple example: 

U " + 7t = O 

which has the solutions u^ ££ 0 and u^ “ sin x on the interval [0,^ ]. 
Both u^ and u^ satisfy the boundary conditions u(0) = u( / 1? ta ') = 0. 

Theorem 6 : Suppose that u^(x) and U£(x) are solutions of (1) 

which satisfy the boundary conditions (2), If h(x) ^ 0 on (a,b), then 

U 1 35 V 

Proof: Set u(x) - u^(x) - u^(x); then u satisfies 

li "+cj(l) 74 ‘ -h h/K) U — o 

with the boundary conditions 

?(fo) - ute) 

By Corollary 2A* u(x) - 0 on (a,b). By analogous reasoning -u satis- 
fies the above equations; hence -u — 0 on (a,b). Thus u S 0 on 
(a,b) . 

In order to incorporate a larger class of boundary value problems 
we consider more general boundary conditions which have (2) as a special 
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case. We consider solutions of (1) which satisfy the boundary condi- 
tions 



y U) COS & -t~ jAn 6> = % 



(3) 



H'(b) Cost + u{b) dm 4 - %L 

where Q- , 4> , % ^ and 2 are P rescri -bed constants with 0 , 

0 i “ . We have (2) when&= 4 ~ *^2 •> 

To derive the general boundary conditions (3) consider the boundary 
conditions 

fi, (~y‘o») +■ 3, v(«) - e, 

At K‘(bj + 3 Z K(b) -- ^ 

where A^, B^, C^, (i * 1,2) are given constants. We may normalize the 



above, i.e., 

Al 



l ?, 



and 



rj Ayr '/ a)) -h ===== 

f At x 4*t ( 4 a?i c? y^+ts, 

A*. . /L) <?, 

+ . /= 7((k) - 

y aa+b* 









71 Ah) 



Then we define the angles & and so that 

Aj. 



Cos 9 ~ 



8, 



K 



and 



Cos <j> - 



y a ?+ xa 

Az 



Cm 6> ~ 



i A 2 + e; 
JL 



Then 0 — and 0 £ <^> ± . 

Theorem 7 : Suppose that u^(x) and ^(x) are solutions of (1) 

which satisfy boundary conditions (3). If h(x) — 0 on (a,b), then 
U 1 ^ U 2 un ^ ess ^ ^ 0, &" » ^ 55 0; in which case u^ and u^ may differ 

by a constant. 

Proof: As before set u = u^ - u^ 9 so that u satisfies 

7f"+ j sx) jv' i-ht/t) 7i ~ o 
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and 

( 4 ) -K'ft) 6&S& -*■ U/&) ^ & - 0 

(5) H*(b) Cos 4 i lifhj $ 1 M <4 ~ 0 

Consider u SL M, a non-zero constant. Then (4) and (5) imply 

M. 9 ^ & - M. <4> = o 

Now M 0; hence sin©- = sin =0 and &■ - - 0. By Corollary 2A, 

M ^ 0 if and only if h *S 0. Thus, if and u^ differ by a constant, 

then h ST 0 and ^ = 0. 

Conversely, assume & - <f> * 0 and h :£ 0; then u SsL M satisfies (4) 

and (5) trivially and 

*yV j/xjw* + fat/) a = o. 

Hence, u^ and u^ differ by a constant. 

Next assume u is a non-constant function which is positive at some 
point in [a,b). By Theorem 2, u attains its positive maximum at a or 
b. Suppose the maximum occurs at x 8 a, Then by Theorem 3, u f (a) ^ 0. 

Now u(a) 0 and 0 £ Q Jfc* , 0 ir <jfi ~ ; hence 

Z "** 

- // '{#) Co JO -h -7/ /nr.) S-u^ & > 0 / 

since cos & and sin cannot vanish simultaneously and are positive 
on the interval under consideration. This inequality is a contra- 
diction of (4). If the maximum occurs at x = b, we have u*(b) ^ 0, 
u(b) 0, and 

?(V6) Cos 4 +~ )> o 

a contradiction of (5). Hence u is a non~constant solution which can 
never be positive. Analogously -u is also a solution which satisfies 
(4) and (5); and ve may conclude that u can never be negative. Thus ti 
SO on [a,b j . 
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Now let the only assumption on the function h(x) be that it is bound- 
ed on the interval under consideration; i.e., we remove the restriction 
that h(x) be non-positive. We will now prove a uniqueness theorem* but 
as the following example shows, the interval considered is of vital im- 
portance. The simple equation 

?i - o 

with u(a) * u (b) = 0 has only the trivial solution for b-a < I'K ; and 
for b-a * uniqueness failed as was shown at the beginning of this 
section. 

Theorem 8 : Suppose that u^(x) and u^Cx) are solutions of (1) 

which satisfy the same boundary conditions (2) with h(x) 0 on [a,b] # 

If b ^ a*, where a* is the conjugate point of a, then u^ SS u 
Proof: Let u ® u^ " u 2 » then u satisfies the equation 

v”+ hit) K —O 

wi Ch 

7//4j =- w/b) - o * 

Since b ^ a* by hypothesis. Corollary 4C implies that there exists a 
function w ^ 0 on [a,b]. Hence by Theorem 4 the function v ® u/w satis- 
fies the maximum principles of Theorems 2 and 3 and their Corollaries f 
Now v(a) = v(b) = 0 since u(a) * u(b) = 0, hence v(x) £ 0 on [ a ,b ] ; 
otherwise, v would attain a positive maximum in (a,b), a contradiction 
to Theorem 2 and Corollary 2A. Since w is a positive function, u(x) £ 0 
on [a,b]. We may apply the same argument to -u since it satisfies the 
same differential equation and boundary conditions, and we conclude that 
-u £ 0. Thus u 2E 0. 
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v . APPROXIMATION IN BOUNDAR Y VA LUE PROBLEMS 



We are concerned with the same boundary value problem as in Section 

IV* Thus, we seek a solution of the differential equation 

(1) (l* ^ ) M - v** hfclu ~ /#/ 

on the interval (a,b) which satisfies the boundary conditions 

( 2 ) 7/mj * ir, ; w (k) - ft. ■ 

Even in many of the simplest problems an explicit expression of 
the solution is impossible to determine. We must then resort to series 
solution techniques or other essentially numerical methods* Obviously 
an approximation solution technique which yields an explicit bound for 
the error compared to the exact solution is desirable* It is possible 
to find such a bound by finding explicit approximation functions which 
provide an upper and lower bound for the actual values of the solution. 

In this section we will develop necessary conditions for the exis- 

tence of such approximation functions, and in Section VI we will in- 
vestigate more closely the approximation functions and related results. 

We assume, as before* that the functions f, g and h are bounded on 
the interval under consideration. With these assumptions it is possible 
to determine bounds for a solution u of (1) and (2) using the maximum 
principle in Theorem 3 without any further knowledge about u* 

Theorem 9 : Suppose that u is a solution of (1) satisfying boundary 

conditions (2) ? and that h(x) 0. If the function Zj(x) satisfies: 

O) ( L +h) l */7 - faj h) 
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(4) 


^ fA) k V, 


u>) k r 7 


and the 


function z^(x) satisfies: 




(5) 


4 J-/ X) 


(»t h) 


(6) 


(H) ir % 4*4 




then 


it M i MU) 


* 


for a ^ 


x £ b. 





Proof: Consider the function 

■VJHO 5 

Then satisfies 

( h^k)[rj3 k o 

on (a,b) with 

V}fo) £ 0 <4 fid VJ(b)&Q. 

Apply Corollary 3A; v^(x) £ 0 on [a,b], Thus, 

MU) * i, /X) . 

Similarly the function 

nri It) - - Mt) 

satisfies 

(li+ft) Cvz J - o 

on (a,b) with 

(a) ± O and fe) 

Thus v 2 ^ x ) — 0 on t a ,b] and we have 

■j??/ k) g. y(Kj — ~£,/x) 

on [a,b]. 

A more general two-point boundary value problem is 

(7) v‘'+ g(t)n'+ kfi) u =-f(t) 

on (a f b) where the solution u satisfies 

- y U) ZOSG + Vrt) s in & * tj 

/ Q \ 

K f (b) <?CS0 + If/h) 2u\4> " t-L 
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The quantities O and are pre-assigned constants wuch that 

06. & 6 on 4 0 6 cf> 6. ^ . 

Continuing to assume h(x) 4. 0, we may extend the proof of Theorem 9 
to get the following. 

Theorem 10 : Suppose u is a solution of (7) satisfying boundary condi- 

tions (8). Suppose also that h(x) 4 0 and & , are in [0,'TV/2] 
and that not all the equalities h 5: 0, © = 0 9 - 0 hold. If the func- 

i 

tion z^(x) satisfies: 

(9) 7 - ~f(£) Oh b) 

(10) (b) (U>scf + -P y (b) S'Uh 4 % 

and the function z^ (x) satisfies: 



(ID 



( 12 ) 



(6+h )[* 2 7 ± 4/t) 

- 'fa) £t)J & 2^ /#) &Pn £ ft/ 

■p/Sb) (U>j 4 + ^b) 4 



then 

t HI*) £r i,OL) 

Proof: Set the function u - z^. Then satisfies 

and 



- % '#0 CoS 0 + Yf fa) ^ & £ O 

fi'fbJcc/tf f vj(h) %^4 -o „ 



If v^ is always positive, then Theorem 2 states that the positive maxi- 
mum of v^ occurs at either a or b. If it occurs at a, then v^(a) ^ 0 
and v^(a) > 0. But 
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- x'(A ) Col 0 + V t (A) 2 -ChQ 6zQ 

can occur only if sin Q - 0, i.e., ^ “0 and v'(a) = 0. Then Theorem 
3 states that v^(x) is identically a positive constant. Thus h(x) 5 0 
by Corollary 2A. Similarly cannot attain a positive maximum at b 
unless 4 = 0 and h 25 0, since 

r; '(b) 60 s 4 + it fb! 3 fa (f 

We may conclude that unless & * 0, ^ * 0 and h S 0, v^(x) £ 0 on 
(a,b); i.e., 

nto) * rfftj. 

By an analogous argument, the function - u satisfies 

{ It+h) [Vi~j ± 0 om 

and 

- Y*) C6>J& -h 7f}(4) &*+ 6^0 

If either h ^ 0 or not both & and are zero, then we may conclude 

(K) ± H/X ) ± -2/fc) 

on [a,b]. 

Remark: If the special case h ^ 0, =0 holds, then from 

the proof of Theorem 10 we have that the functions and are both 
identically positive constants, say 

& M ,> o and v± U)= J% z > o. 

Then 

-2,(1) <2,Wi-M, - UU) - £>/*)- 

i.e., we have 

2 , Uj < y/x) < 2} ( x ) 

and z^ and z ^ differ by a positive constant. 
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We drop the restrictions that h 0 and that & , ^ are in [0, / 2 ] • 

From the discussion on the background of the generalized Sturm-Liouvil le 
boundary conditions we may have just as easily taken €> , ^ in (- / ft*' /2/ftV2] 
with no loss of generality. 

Theorem 11 : Suppose that u is a solution of (7) satisfying boundary 

condition (8) with - ^ /2 ^ — ^72, - / ^/2 ^ /2. Let 

the functions z^ and z^ satisfy the inequalities (9), (10) and (11), (12) 
respectively as in Theorem 10. Then the bounds 

y(tj & 

hold in (a,b) if and only if there exists a positive function w on [a,b] 
satisfying 

<i3) ( JLi't'k) [nj~] — 0 fat b) 

- 14/ ) <*&j q + bs*. & ^ c 

(14) Uj ( (h) q#j cvfb) Sih<f ^ 0 

in such a way that not all the inequalities in (13) and (14) are equalities 
Proof: Suppose there exists a positive function w on [a,b] which 

satisfies the inequalities (13) and (14). Then define the function v = 
u/w. As in Theorem 4 

(W)[tO = ( L-H,) 

- v"-b[z f-' + = £ 

which we may rewrite as 

(is) ( j A -hrt) [Y] - r v + tfr' + M ir j-j 

where G = (2w“/w) + g and H - l/w(L-*-h) w. Note that H, the coefficient 
of v, is non-positive. The functions v satisfies the boundary condi- 
tions : 
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-j yfa) mj +rf*j wtoj] + vto) m&) 
jy f (b) Uj/b) t v(b) u/'ftyjdc \S<f +■ 11(b) Uj{bJ #wi 4 ~ 'ft?. 



which we may rearrange to form 



(16) 



-y^JWIW e f i/fflf-w'tyQeie-i w/M) ©] « T, 

vh)»/(b) ■*■ tfk)hv'(h)aoi4± t v(y)&K<f] • ?i 



Dividing through by w(a)cos <£?■ and w(b)cos^» , then multiplying by the 
non-negative quantities cos <5* and cos ?* in the first and second condi- 
tions of (16) respectively, we have 



-vh)Ccj£ + M) 



COdG (-to' fa) 

U/fGd Cgjq- 



CoJG 



(17) 



r'ft) <%sj> + i^b) 



ecu? ( to '(b) Cos </ +U/f bj 
U)(b) Coicf* 



4 Z (hs<fi 
<v(b)eGS<fi 



where the angles & and ^ are defined by 



© s. 



Case (~ wYa )<%>$& + u/ fa) ■& ) 

W4) COG 






Cos 4> (u/'ft) ■+ Uf(b) Hk&J 

b//b) CoS ft 
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or equivalently 



-/a 



H O 



— Ui'tA) COS& -f tv&) Shi & 

(v (a ) Cosq 



u/( t*- 

7v(b) Cos 4 

By the hypotheses on & , 4 and w, both tan d> and tan <fi> 
negative. By (14) we may always take <j5 and jf' such that 



(18) 



7 ^<W]<4> 



are non- 



Ot €> < ' f y z 0±.<f - - 



Note: If 9- = >72, take = IV /2 ; if <^ = 7^/2 take £ = ib / 2. 

Hence we have 

-r'htjcose -t i/fa) 

(19) ___ 

X ( (b) cos 4 ** vOi) $<h*<p = T x 

with © ,4 in [0, TV /2] and ^ ^ cos /w(a)cos© and ^ = 

cos ) /w(b) cos 

If the functions z^ and z^ satisfy the conditions (9), (10) and 
(11), (12) respectively; then the functions z^/w and z^/w satisfy the 
analogous conditions with respect to equation (15). Hence by Theorem 10 
the following inequality holds 

2 % U) ^ W/l ^ *,(* > 

tis(Z) " tO(Y) 0U(l) 

unless 0> = ? ~ 0 and H S 0. From the proof of Theorem 10 this situa- 
tion holds when the inequalities (9), (10) and (11), (12) are equations. 
Hence, if there is a positive function w(x) satisfying (13) and (14) but 
such that not all are equations, we may conclude 



( 20 ) -2Jy) £ 7i(X) 
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Note: If w satisfies the boundary value problem 



(L+h )M ~ 0 ( tit b) 

! - HJ'te) -t uj/cu 3wi6 « O 

lAj’(b) a>$$ -h w(b) 2^4 ~ 0 

then let ^ be any real number and u a solution of (7) which satisfies 
(8). Then 

(li+k)[ni t§v]« (MJltf + fU+WM =o 



and 




+ 





&-M & — 



-b [ u( b) b f 







Thus, if w satisfies (13) and (14) as equations, then any multiple of 
w, when added to a solution of (7) and (8), will yield another distinct 
solution. Hence, if at least one solution u exists, then there are 
many. 

Now consider the boundary value problem 

(2i) (It't'h) l.w] “ 0 Pyi (#ih) 

with 

i - iVY 4 )C£>S& wfa) - 1 

Co i/b ) Cos 4 + cvfo) 3 m 4 ~ l 
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By applying the classical existence theorems, let w be a solution of (21) 
satisfying (22) . Noting that the function w satisfies analogous condi- 
tions to the function u in Theorem 10, we have 

± W(%) £ - 2 , (KJ 

for appropriately chosen bounding functions z^(x) and z^(x)* Also 
z^(x) 0 satisfies condition (11), (12); hence the function w(x) is 

non-negative on [a,b]. We now show that w(x) is actually positive on 
[a,b] . 

If w is zero at an interior point, then w* is also zero there, 
otherwise w would have some negative values* But the identically zero 
function satisfies the same initial conditions as w(x) and by Theorem 
5 for uniqueness of initial value problems, w 5E 0 on [a,b] which con- 
tradicts (22)* Hence w cannot vanish at an interior point. 

If w vanishes at an endpoint, say at a, then from the first condi- 
tion of (22) we have 

(aj Y*) ccs & “ - i . 

Now w non-negative implies that w'(a) £ 0, a contradiction since cos<£*- 
^ 0. Thus w(a) 0. Similarly w(b) ^ 0, and we have w ^ 0 

on [a,b]. 

As with Theorem 1 and Corollary 1, it is often possible to prove a 
more general theorem and by a judicious choice of parameters obtain a re- 
lated, more specific result. For example, if h - 0, 0 ^ /2, 

0 £ £ *1^/2; then the function w zz. 1 satisfies (13) and (14) and 

Theorem 10 follows immediately from Theorem 11. 

The fact that the function w does not appear in the inequality 

€ Z CKJ 4? 7fte) £ fx) 
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prompts us to try to obtain a necessary and sufficient condition that 
guarantees that the functions z^ and form upper and lower bounds for 

u independent of a function w. 

Theorem 12 : Suppose that the functions z^(x) and satisfy in- 

equalities (9), (10) and (11), (12) respectively, in such a way that 
equality does not hold in all the conditions. Let g(x) be bounded be- 
low on each interval [a,c] and bounded above on each interval [c,b] 
with a ^ c ^ b. Let u(x) be a solution of (7), (8). Then the bounds 

(23) ^/K) * 

hold if any only if z^Cx) ^ z^(x) for x ^ b. 

Proof: Trivially, if (23) holds, then z^(x) ^ z^(x). Hence assume 

~ Z 2 non ” ne g at ive and we must show that (23) holds. Put 

= -z ( frj - zjt) 

If q ^ 0 on [a,b], then q satisfies the conditions for the positive 
function w in Theorem 11; i.e., 

(24) (L4hJ[% 7 = Oh fy,b) 

since z^ and z^ satisfy (9) and (11) respectively. 

- [-?/«) - Cose + - '&4)j ko 

(25) fe'/hl ~ Cos 4 [$/*)-*/*}) ^0 

since z^ and z^ satisfy (10) and (12) respectively. Hence we may con- 
clude by Theorem 11 that 

■2JX) £ Vff) £ - 2 , (>C i . 

Thus, we need only consider the case where q has a zero in [a,b]. We 
will show that q cannot vanish at an interior point and that if q 
vanishes at one or both endpoints then we may find suitable positive 
functions w(x) such that Theorem 11 will hold. 
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Suppose that q(c) = 0 for c in (a,b). Then q'(c) = 0 since q ^ 0 
on [a,b], i.e,, q has a local minimum at c. We may conclude that q ££ 0 
on [a,b] by the uniqueness of solutions to the initial value problem. 

But then equality holds in (24) and (25), a contradiction. 

We are left with q ^ 0 interior to [a,b] and q vanishing at an end- 
point, say at x - a. Then q*(a) ^ 0 by above. Thus ££ » / fK /2, other- 
wise the first condition of (25) would be violated. Similarly if q(b) * 0, 
then q f (b) 0 and <^> - /2. 

We now consider the case where q(b) - 0 and q(a) ^ 0, and proceed 
to find a positive function w such that Theorem 11 will hold. From the 
preceeding paragraph <& = ip /2. Consider the initial value problem 



(26) 



(u+h)bti = 0 

■}'fy)~Ccs& ansi 



which has a solution by classical existence Theorems for differential 
equations. Let r be the desired (unique) solution. We note that for 
B < /2, r is (strictly) positive at a, and if 0 * 'V /2, r(a) 

* 0 and r is positive near a. 

Now define the function v = r/q for a < x < b. Also note that 

rfo) v 

nrm ~ — ^ o 



and 



Jfa) 8 - jp'fW e&i© 

J>#j 



by the first condition of (25). Using the facts that (1+h) ( r] = 0 
and (L+h) (q) 0, we may conclude that v satisfies 



= V ,! i NT = 0 
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with G = (2qVq) + g ^nd H » l/q(L+h) [q] ^ 0, as in the proof of 
Theorem 11. The coefficients G and H satisfy the restrictions in Theorems 
2 and 3 on any subinterval of the form [a,c]. Hence either: (i) v(c) 
y v(a) and v' (c) 0, or (ii) v(x) S v(a) for x in [a,c]. In either 

case r(x) y 0 for x in [a,b) since both v and q are positive. 

If v(x) is identically the positive constant r(a) on [a,b] , then q 
is proportional to r and hence must satisfy the same equations. Thus 

7 = O 

and 

-yfa) C$3 & + r/a) s^yi G ~ ~5'm$Oose +sine Ccso - 0 

implies that 

- ease -h ffa) & - o. 

By hypothesis on z^ and z^ (hence q) equality does not hold in all 
conditions (24) and (25). We must then have 

3*(b) c&i -h 2(b) > O 

s' * / 

from the second condition of (25) and that 1 2. Hence q(b) 0, 

a contradiction. Thus we are left with case (i) , i.e., for some number 
c in (a,b), 

vte) > ~v(a] ^ o v'(K) q -Jcr ( . 



Now 






r(X) = 

Put ^ (x) = r'q - rq'. Then from (24) and (26) 

(2?) - r ~° 



and 

( 28 ) f = = <a 
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Solving for -qhr in (28) and substituting into (27) we see that 
satisfies the differential inequality ^ 7 — 0. 

On the interval [c,b] g £ M by hypothesis, thus 

Solving for fi we may write 

_ JA(x-c) 

fie 

and for x - c 

He) ± A J 

thus 

-Mfr-c) 

fit) i fie) e 

In particular 

_ M (X-e) 

fi(b) > ft(c) e > 0 

since v(c) 0. i.e., ^ 

r'tkljj.lb') - j'(h)r^b) Ao . 

But q(b) = 0 and q’(b) £r 0, therefore r(b) 0. Hence we have r 

positive on (a,b] and satisfying 

U - O 

with 

-r'w cos (9 * r/aj sw e 

and since <£> - ^ /2 

7*Y6) <fi = yY b) > c . 

Hence the function 

w/x) s yt) + v(k) 
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satisfies the requirements oi Theorem 11, and we have 



U(X) - 

On the other hand, if q(a) - 0 and q(b) 0, analogously we 
may show that the solution of 

C = 0 

l ;0>) - Cts4 and sY b! - 

is positive on [a,b), so that the function 

ti'/XJ s fU) V 3( %) 

satisfies Theorem 11. 

Finally if q(a) = q(b) = 0, we have r ) 0 on (a,b] and s 0 
on [a,b) with r(a) = s(b) = 0; hence 

rffl+stx) 

satisfies Theorem 11. This concludes the proof of Theorem 12. 
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VI. APPLICATIONS OF APPROXIMATION FUNCTIONS 



Through the preceeding sections we have seen that the results have 
been essentially theoretical and that examples have not always been 
fruitful or plentiful. Sections 1-4 present one -dimensional maximum 
principles and their use in some differential equations theory while 
Section 5 actually gives some practical techniques to apply in boundary 
value problems. The problems we are considering are second-order dif- 
ferential equations with variable coefficients which in general must be 
solved by series or approximation methods. 

There does not seem to be a ’"best" form for an approximation func- 
tion since, in general, each boundary value problem needs individual at- 
tention. However, there is a general procedure to follow to determine 
bounds and behavior patterns of the solution (under certain conditions) 
and then the approximate value of the solution at any point in the inter- 
val under consideration: 

1. Find the straight line z(x) which satisfies the boundary 
conditions as equalities and apply the operator (L+h) to it. If (L4h) 

[z] is of one sign on the interval, then z will serve as either an upper 
or lower bound as determined by the sign of (L+h) [z]. (i.e., (L+h) [z] 

0 implies that u(x) £ z(x) and (L+h) [z] >0 implies that z(x) 
u(x)) . 

2. If the constants ^ ^ Q f the boundary conditions are 

of the same sign, then the approximation function z(x) Sr 0 will show 
that the solution u(x) is either non-negative or non-positive. (i.e., 

^ 0, ^ 2 — 0 implies u ^ 0 and ^ 0, ^ 0 

implies u ^ 0) . 
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3. If the general behavior of u(x) is desired, it is only 
necessary to bound u between the straight line from (1.) (assuming this 
is possible) and some constant function which satisfies the approxima- 
tion boundary conditions (9) and (10), or (11) and (12) in Section V. 

4. If an approximate value of the solution u(x) is desired 

at some point in the interval then find a polynomial of exponential func- 
tion that satisfies the approximation boundary conditions mentioned in 
(3.) above. If a more exact approximation is desired, then the form of 
the approximation function is more complicated and the difficulty of evalu- 
ating and checking the restrictions increases. In many instances it 
would be as easy to find an approximate or series solution (if initial 
conditions are known) and evaluate it numerically on a computer. The 
value of the above technique is that the behavior of the solution of a 
complicated second order differential equation can be determined easily 
and in less time than other methods. 

The following examples will help point out the applications of 
Section 5 and Theorem 12 specifically. The choice of the approximation 
functions may seem ambiguous, but they readily become more understandable 
in the context of each problem. 

Example 1 : Consider the boundary value problem 

(1) ?/' v / 3V- XU ~ Q OH ( '0,0 

wi th 

( 2 ) - 21 (t>) “ 0 tfftd 7/ 0)~ t . 

As noted in (2.) of the preceeding discussion the solution of (1) satis- 
fying (2) is non-negative. Consider the straight line ~ ^ (Ki'l) 
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which satisfies (2) * Then 

= i [3~ £ 0*xf] ± 0 

on (0,1). Hence let z(x) = z^C*) and we have 

{(%*!) i VOC) , 

We note that y 1 satisfies the requirements for an upper bound, i.e., 

O 

(L+h) [y] ^ 0 and y(l) = 1, -y* \0) + y(0) = 1 <> 0. We have quite 
simply determined the general behavior of u in (0,1), i.e., 

(1/2) (x+1) u(x) £ lo 




Fig* 2 

To determine more accurately the bounds for u at any point x, we choose 
a more complicated function for z^, the upper bound* A natural choice 
to consider would be 



z j(x) * (1/4) (x+3) 

which satisfies 

z(l) - 1 and ~z ? (0) + z(0) = 3/4 > 0* 
But z(x) does not satisfy 

o> (M)i3, 7 - i[i- *'*+*)] 
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on the entire Interval (0,1), Bat we may define 



{g[& 4S lj 0ix *- i/vi-s] 

%n) = J 

\i (f*') 3 i[^\-s]cx&i 

where (1/2) ( V 21-3) is a zero of 3 - x(x+3) = 0. Then z^ is continu- 
ous and satisfies the conditions for an upper bound and we have 

^ (K+i) £ Ml) £ £ [aTT *bJ 

on (0 , (1/2) ( *T?l - 3)) and 

~2, (t+i) ± h(k) - ifiHS) 

on ((1/2) ( fn - 3),1). 

Example 2 : Consider 

(4) u M - (l+x^)u = 0 for 0 < x 1 

with 

(5) u(0) = 1 and u(l) « 0. 

The solution u must be non-negative. Consider the straight line 

z 1 (x) = 1-x 

which satisfies (5), Now 

~ — (H&*) 2/ 

* 0 — <7yx 2 ,) rV-x) ^ O 

since 0 <£ x 1. We now have an approximate idea of the shape of u. 

We pick for (the lower bound) a parabolic function. We choose z^ 

such that * 1 and z^l) = 0, and so that is concave upward. 

In Example 1 we showed steps 1-3 of the general method and here we show 
how a better bound can be found by choosing a more complicated approxi- 
mation function. Let 

z 2 (x) * (x-1) 2 
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which satisfies (5). Then 



Hence 



- (/*x*)*x 

- Z- - ° 0,1 foO ■ 

(t-i) i £ VK) £ /-X 



In particular x - 1/2 

j £ H({) £ 

Choosing for the expotential function 

€*(>-*) 

which satisfies 

%.(<>)= I and 




Now 



(Irfh) [$i J~ ■?2 / — 



($-*) <t K - (/+L l ) O- X)* * 



* £ 



on fd,t) 



-x 3 2 

since e 0 and 2 x - x „ Hence 

i* a-*.) ^ «<*■) *= /-x, 

and in particular at x * 1/2 



.303 £ u(x) £ .500, 

a considerable improvement. 
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Since z^(0) ~ z^CG) ® u ^)» we may bound u^O). Now 

*'m - ^ f^r^! 

ic~> o i * J ) 

and from our previous results 

-2 - Hi / o ) ^ H f fo) - to) * ~*1 * 

c 

Example 3 ; Consider the problem 



( 6 ) 



V* ~X7i -O Ttn 0 <X < / 



with 



(7) 



- v'(c) -t m/c) and ~u(Q ~ / 



Noting the similarity of boundary conditions (7) with boundary condi- 
tions (2) in Example 1, we choose 

?/*)=. i (K-tt) 

and note that 

a**)[gj- 9 ¥ ~*i 

= -4 ttt+O-o w (iA). 



Hence z^(x) ~ (1/2) (x+1) is an upper bound for the solution u(x). We 
choose for z^ the exponential function 



e 

which satisfies (7). Then 



*W 






= 0-t>c ^ o t>*> fa,t) . 



Thus 



£ HU) 4z. ■£ (*.+/) j 
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and in particular at x = 1/2* 

0.6065 ^ u (1/2) ^ 0.7500. 

An interesting application of the theory developed allows us to 
bound solutions to the differential equation 

(jj+h)M “ h(K)n -■ iw 

on (a,b) even when the coefficients g(x) and h(x) are piece-wise continu- 
ous. 



Example 4 : Consider 

( 8 > hnc)K — o ~fcy < 

with 

( 9 ) y/oj-Q and O) v / 

where 



1 ; 04x il/2 



h(x) = 



2 ; 1/2 < x & 1 
From the boundary conditions consider the approximation function 



z(x) * X 

which satisfies (9). Then 

hr*) 2 - 



K-t -O y 



Hence 0 £ u(x) x» since u must be non-negative by (9). A better 
lower bound may be determined by considering the parabola 

, v 2 

z 2 (x) = X 

which satisfies (9) and 
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= if + - hot) * 



X l + 2 >O r ) O'ZXSr *i 

2 J 



Thus 

x 2 £ u(x) £ x. 

When working through the theoretical background for this technique , 
there is a question of generalizing the approximation functions so that 
the "best" bound for the solution u can be obtained. In considering this, 
the following discussion will point out some problems which arise. 

We will choose general expressions containing parameters for the 
approximating functions z^ and z^ and attempt to pick values of the 
parameters to make z^ and Z 2 satisfy the appropriate conditions* For 
example choose _ ^ ) n 

Cl) ~ ft I 2 ~ € J 



and 



4 ) = 8 [z 



e 



1 



where A,B and &( are constants to be determined. From Section V we 
know that z^(x) must satisfy 

(7^J/2/7 ~ on (a ,b) 

i, (4) h7, ) 3 (b) 2 £ ■ 

( L ^h) [2/ 7 ~ ° + *} fXJ'2) 1 'Z) 

sr - A £ C< fX 



and 



Then 
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Now choose so large that 

“ + hit) > Oj 

where f,g and h are assumed to be bounded on [a,b]. To make z^(x) an 
upper bounds we must have 

(Mh)litl ~ W htt) - A (< - ft) -hh/tt ) e J ± ^/X) ) 



or 



4 1 2 hit) - ( os 1 - *< «j {K) + h M) e J ^-fvc) 



This is satisfied if we let 



k * /mm /- 2^(1) + A /2 ~ uq tx) 4 6//j ) e 

jK A \ r Z. i {-S 



ot/X-4) 



di K&k> L 

and then if we choose A such that 



fi ' j: j- ~h*>\ . 

By its definition k ^ 0. We have now satisfied the condition 

(M)fe 1 A^ftx) acy < t>. 



We must also satisfy 

a 



• A) 



-({$)- A j 2- <? J ~ A ^ ?7 j 

r — •“* ) 7 v. — 

*/ h) ~ A £ 2 - tf 7 — . 



Let c - 2 • e > and choose A such that 

A- ) — ) i *«*■}-■(*)} ) O] 

i.e., A is the largest of the four numbers in the brackets. 
Then 

u(x) z 1 (x). 
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Using the same in z^(x) , we choose B as the smallest of the four 
numbers 



B * min 



\% , £ 5 — V v o7 

L c. ' ^vcfri ^ J J 



; 



then 



%[z-e J & A[2-e J 



on [a,b]. 

Example 5 : We solve the following problem in two ways. Consider 



with 



u" - xu = 0 for 0 < x < 1 



u(0) = 0 and u(l) = 1. 



Let Zj(x) = x which satisfies the boundary conditions, and 



(L+h) [z^J = -x £ 0 on (0,1). 



Hence 



0 u(x) ££ x. 



For z ^ we choose 



=s 

which satisfies the boundary conditions, and B ^ 0. Then 

(h-Lhjfet]- 0 on (0,1) 

for B ^ 1/2, which can be seen by considering the values of B which 

make 

o 




since 



/ 



(M) HO I > 0 , 

U-0 
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and we are dealing with a cubic equation* Hence we choose B = 1/2, and 
we have 



(l/2)(x-fx^) £ u(x) x^ 

since z^ ib z^ on (0,1)* From the theory developed before this problem, 
choose the approximation functions as follows: 

A[z-e^] 

and 

&[Z- 

Then we must choose so large that 

% >C on (0,1), 

hence choose CX - L Then 

k = min ^ = 1 (at x = 0) , 

and 

9 - 1 

c * 2 - e 

Hence 

A = Max \ 0) \ 0 j (0 I cr- 

L z~e" 1 J 2e.-\ 

and 

B = min * O- 

Hence 

o £ mi.) & — £— / ^ ^ 

which does not improve on our previous results* 

The method outlined previously and illustrated in the example also 
leads to very interesting aspects of bounding solutions of second-order 




57 



boundary value problems with general boundary conditions. 



We may find the straight line 

do 2/*;= 



which satisfies the general boundary conditions: 

/- n‘{a)eos & + 2(4) s<h & * ^ 

(ii) < 

(. + 2(b) 5iVi 4 “ ^ 

if the parameters $ and satisfy certain restrictions. 

We may rewrite (11) as 

a, (- nos e) -a «* - 37 

c*» <?£>i ^ v- (0, b + <?») 

or equivalently 

C t l# & - aosoj + sin & - 

+ oos <f] + ‘ft* ; 

which has a non-trival solution for c^ and c ^ if the determinant of co- 



efficients 

Z> = 



-cose) 




(b 5 

is non-zero (assuming "^l 2 * + 0) • 

x> s 5/« a 5/> ^ - [*•'” 9 Cos/* + CG&Sin/J 

- {4-b) Sin 6? — s/>? (o+&) . 

Hence if and ^ are not solutions of 



(12) *sfw(& + /) — s/>ie 5 ^t / 
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n'TV 



then we may always find the line (10) * We note that if = 

and O = nl^ 9 <$> “ rilV for (n * 0, 1,...), then (12) is 

satisfied. 

In Example 1 we used a polygonal line as an upper bound. This 
procedure may not be fruitful in the non-homogeneous case. The reason 
for this will be made clear later. 

Example 6 : Consider the non~homogeneous boundary value problem 

2- 



with 



'll" - 71 - Zit-O 
7((0) - 14{/) « o 



To determine possible concavity of the solution, consider 

ZtK) = 1 

then 

* -/ < 2Y1C-X 2 ’ CH fcjO) 

hence z = 1 is an upper bound and the solution u(x) is concave upward 
since 

2" = - I 4 2 <0 (Oji). 

To determine a lower bound, consider 

2/X) - 3 

then 

(M)[X- 2-2 = -5 >•?-'*-**' On (0,1)', 

hence z » -3 is a lower bound since the appropriate endpoint in- 

equalities hold. 

By substitution we may verify that 

lOxJ - p 

is the solution. We might expect then that the polygonal line 
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- x ; 

— j fy^X * 
i-i j ?<f < x ^ i 







would satisfy the theoretical conditions for a lower bound since 

z^(x) — u (x) on (0,1). 

0 ^ % 1 




Fig, 4 

Note that this is the opposite of the procedure in our theory where we 
determine bounding functions by applying the restrictions to them 
rather than determining whether actual bounds of the solution satisfy 
the restrictions. Recall that for z to be a lower bound* then 

(i) (L-fh) [z^] — f(x) s 2 + x * on (0,1) 

(ii) z^CO) — 0 and z 2^ “ 

Applying the operator (L4h) to z^ 9 

L i 

^z-nc-x 3 ' 6* (o,0 , 
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but 



^2 Jo) - o a,v\ d 2ji) = 0 



From our theory this would seem to imply that z^ is an upper bound for 
u; but by an oversight we have not considered the fact that each segment 
of the polygonal line must satisfy the restrictions individually. 

Hence for z (x) = -x [0 9 l/4] we see that 

with 

2jc) t&vtJ *±(ij = 

Thus our theory does not apply. We find similar results for the re- 
maining segments of z^. 

This points out an important result that for the non -homogeneous 
case the theory developed may not apply to the polygonal line approxi- 
mation functions. 

The reason for this is clear upon consideration of Example 1, where 
we choose the straight line 

iW 1 i (l+3) 

as an upper bound s and find that 

(l, * h)l 2 2 3* - K £ = y [$ “ ^ C 

changes sign on (0 9 1), i,6o, (L+h) [z] has a zero in (0 9 1) and by 
suitably changing z(x) we can satisfy the restrictions 

(Li+h) - O <9H (Ojf). 

But for the non-homogeneous case 9 we must satisfy 

(/j + ~f(K) 

and we may not be able to suitably change our approximation function 
to satisfy the conditions. 

This points out the importance of parametric approximation func- 
tions and being able to choose the parameters in such a manner as to 
satisfy the restrictions. 
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Consider then 



3 nj~ ft [ 2 - e *3 
ij/W * 6 \z - e*'] 



and choose A, B and o( according to the theory developed preceeding 

Example 5 so that z^ and z^ are upper and lower bounds, respectively. 

2 

Then must satisfy ^ « 1 ) 0; hence, choose ^ = 2. Then 

choose 



k** 

t 



n Tl + $€ J = 2iS£ 2 X ^ £. 



We must now determine 



mi S-'ftn f - jx*-x~ z j 

0£f*l ‘ ' 0-x.£. ( 



~ -2 rrK=)'j 



and 



Him i~ 




~ flt-VK* 



* - a ^ ^ ' 



Since ^ 35 0, we do not need c ~ 2 - e ^ since ^ 9 /c s 0. Thus 



choose A and B so that 



A = nwi / £, 



-2 






2 ^ 3s? 



- Z 



>0=0 



and 



3= , 15 

1 i-mz-i- > 



2+ a e 



'Z 



Thus 






£+ /2. £ 



= a~"j 



a 'y(t) £ o „ 
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We may graph the above by approximating the function z^(x), and we 




Fi go 5 

We now have the option of attempting to change the parameters A, 

B and ^ to obtain better bounding functions and or keeping the 
existing functions and tolerating the possible large error* In either 
case we are sacrificing time or exactness and are defeating our pur- 
pose,. 

The technique used in the proceeding examples has indicated the 
possibility of approximating solutions u(x) by multiplying the straight 
line s(x) found above by ^warping 5 ” functions ^ (x) and ^ (x) so that 

£ K(t) - 5U) 

This technique should not only help determine the general behavior of 
the solution but yield usable numerical approximations* 
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